Electronic properties and magnetism of iron at the Earth's inner core conditions 
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We employ state-of-the-art ab initio simulations within the dynamical mean-field theory to study 
three likely phases of iron (hexogonal close-packed, hep, face centered cubic, fee, and body centered 
cubic, bee) at the Earth's core conditions. We demonstrate that the correction to the electronic free 
energy due to correlations can be significant for the relative stability of the phases. The strongest 
effect is observed in bcc Fe, which shows a non-Fermi liquid behaviour, and where a Curie- Weiss 
behaviour of the uniform susceptbility hints at a local magnetic moment still existing at 5800 K 
and 300 GPa. We predict that all three structures have sufficiently high magnetic susceptibility to 
stabilize the geodynamo. 
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PACS numbers: 

Being the main component of the Earth's core, iron 
attracts significant attention from broad research com- 
munity. Understanding its properties at ultra-high pres- 
sure and temperature, ranging from studies of the core 
structure to modeling the geodynamo is a long-term goal 
for the condensed matter physics, and is essential for ex- 
planing geochemical observations, seismic data, and for 
the theory of geomagnetism, to mention few examples. 
In spite of all previous theoretical and experimental ef- 
forts, the crystal structure and properties of sohd Fe at 
the inner Earth core conditions remain a subject of in- 
tense debates. All three phases stable at low pressure- 
temperature conditions, namely, hep, fee, and bcc, have 
been suggested as possible crystal structures of iron or 
its alloys in the Earth inner core [ll-@|- 

Theoretical simulations of iron at high pressures and 
temperatures generally rely on the picture of a wide-band 
metal with insignificant local correlations Indeed 
under compression the overlap between localized states 
increases and so does the bandwidth W, while the local 
Coulomb repulsion U between those states is screened 
more efficiently. The reduction of the U/W ratio is used 
to rationalize the absence of electronic correlations be- 
yond the standard local density approximation (LDA) at 
high-pressure conditions. The increase of the 3d-band 
width also results in the corresponding drop of the den- 
sity of states at the Fermi energy leading to disappear- 
ance of the driving force for magnetism according to 
the Stoner criterion. Besides, even at the ambient pres- 
sure but at very high temperature T » Tc (Tc is the 
Curie temperature, which is 1043 K in a— Fe) local mag- 
netic moments are expected to be suppressed due to one- 
electron Stoner-type excitations. Thus, when extremely 
high pressure and temperature are simultaneously ap- 
plied, disappearance of the local magnetic moment seems 



to be inevitable. Due to these considerations iron at the 
Earth's inner core conditions has been modeled as non- 
magnetic within LDA-based approaches. The results of 
recent works of Sola et al. 0,01, who applied fixed-node- 
approximation quantum Monte Carlo techniques to com- 
pute the equation of state and the melting temperature 
of hep Fe at extreme conditions, are in good agreement 
with previous LDA-based simulations, thus strengthen 
the above argument, at least in the case of the hep phase. 
On the other hand, Glazyrin et al. Q have just demon- 
strated importance of correlation effects in hep iron re- 
vealed by an electronic topological transition induced at 
pressure 40 GPa and room temperature. 

So, are electronic correlations important at extreme 
conditions, and can they lead to qualitatively new phe- 
nomena? To address this question we have investigated 
the impact of correlations on the electronic structure, 
magnetic properties and thermodynamic stability of iron 
by performing ab initio simulations of the bee, fee and 
hep phases for the volume of 7.05 A/atom, correspond- 
ing to the pressures expected in the inner Earth core, and 
temperatures up to 5800 K (The c/a ratio in hep Fe was 
fixed at 1.60 [lO|). We ernploy a state-of-the-art fully self- 
consistent technique [111, [Si combining the full-potential 
linearized au gme nted plain-wave (FLAPW) band struc- 
ture method [13| with the dynamical mean-field theory 
(DMFT) [Hi treatment of the on-site Coulomb repul- 
sion between Fe 3d states. A combination of LDA and 
DMFT has been applied earlier to study thermodynamic 
stability [l^ and to describe magnetic properties [l^ of 
paramagnetic bee Fe at ambient pressure. 

In Fig. [1] we display the LDA-I-DMFT k- resolved spec- 
tral functions ^(k, E) for the three phases obtained for 
the temperature of 5800 K. First, one may notice that 
in hep Fe the electronic states in the vicinity of Ep are 
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FIG. 1: (Color online) The LDA+DMFT k-resolved spectral 
function A(k,E) (Vat/eV) for bcc (a), fee (b), and hep (c) Fe 
at volume Vat =7.05 A'^/atom and temperature 5800 K. A 
non-quasiparticle Eg band is seen in the vicinity of the Fermy 
energy along the N — T — P path in (a). 
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FIG. 2: (Color online) The ratio of the inverse quasiparticle 
lifetime F to temperature T vs. T. The solid red, dashed blue 
and dash-dotted green curves correspond to 3d states in fee, 
bee, and hep Fe, respectively. They are split by the crystal 
field into t2g (diamonds) and Sg (circles) representations in 
the cubic {bee and fee) phases, and two doubly-degenerate 
{E' and E" , shown by diamonds and squares, respectively) 
and one singlet {A'l, circles) representations in the hep phase, 
respectively . A non-linear behavior of F/T for bee Fe eg 
states is clearly seen. 



sharp (their red color indicating high value of A{\i,E)), 
hence e— Fe exhibits a typical behaviour of a Fermi-liquid 
(FL) with large quasi-particle life-times in the vicinity of 
Ep. In contrast, the bcc phase features a low-energy Cg 
band along the N — T — P path that is strongly broad- 
ened, thus indicating destruction of quasiparticle states. 
fee Fe is in an intermidiate state, with some broadening 



as Eb = D{E){E - C)dE, where D{E) is the corre- 
sponding LDA partial density of states (PDOS), C is the 
ccntralweight of the band. Resulting E^ for the Cg and 
t2g bands in the bcc [fee) phases are equal to -1.05 (-1-01) 
and -1.08 (-1.20) eV, respectively. One may see that the 
difference in kinetic energy between the Cg and t2g bands 
in bcc Fe is rather small and, in fact, even smaller than the 
corresponding difference in the fee phase. Hence, it can 
hardly explain the observed qualitatevely distinct non- 



noticable in the Cg bands at Ep in the vicinity of the F FL behaviour of the Cg states in bcc. It has been pointed 



and W points. 

To quantify the degree of non-Fermi-liquid (non-FL) 
behaviour we have evaluated the inverse quasiparticle 
life-time F = — Z5 [I](iO+)], where the quasiparticle 



residue Z 



= 1 
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-Q^ — itj^-u+j by extrapolating the 
imaginary-frequency self-energy Yj{iu)) to zero. In the 
FL regime F scales as T^, hence F/T vs. T is linear. In 
Fig. [2] we display the temperature evolution of F/T for 
the relevant irreducible representations of the Fe 3d shell 
in all three phases. One may see that F/T in hep Fe 
exhibits a linear increase typical for a FL up to tempera- 
tures expected in the inner Earth core. In contrast, F/T 
for the bcc iron Cg states features a linear and steep rise 
for T < 1000 K and then behaves non-linearly, indicating 
a non-coherent nature of those states at high tempera- 
tures. The bcc Fe t2g and fee Fe Cg electrons are in an 
intermediate situation with some noticeable deviations 
from the FL behaviour. 

The tendecy of bcc Cg states to a non-FL behaviour has 
been noted before for ambient conditions and explained 
by a smaller effective bandwidth of the "localized" Cg 



band compared to the t2g one [17|. We have evaluated 



the one-electron kinetic energy of the Cg and t2g bands 



out [18|, [19| that a van Hove singularity in one of bcc Fe 
Cg bands leads to formation of a narrow peak in the cor- 
responding PDOS in the vicinity of Ep (see SI Fig. SI). 
A large peak in PDOS located at Ep leads to suppression 
of the low-energy hopping and to the corresponding en- 
hancement of correlations, as has been recently pointed 
out for the case of Sr2Ru04 [23|. A similar suppression 
is observed in the Cg hybridization function in bcc Fe (see 
SI Fig. S2). 

Having demonstrated the impact of correlation effects 
on the Fe electronic structure we now focus on its conse- 
quences for the Fc phase stability and magnetism. To 
evaluate the impact of correlation effects on the rel- 
ative stability of the three phases we have computed 
the corresponding correction to the fixed-lattice free en- 
ergy by employing a coupling-constant integration ap- 
proach (see in, e.g., Ref. [2l|). We define the free energy 
F\ = — i InTr (exp[— /3(iJo + A77i)]) corresponding to a 
given value of the coupling A G [0 : 1], where Hq is the 
one-electron (LDA) part of the Hamiltionian, Hi is the 
interacting part equal to the difference between the Hub- 
bard term Hu and the double-counting correction E^c- 
The coupling constant integration leads to the following 



expression for the many-body correction: 

AF = F-Fo= [' ^-^^dX (1) 

In derivation of Eq. [1] we neglected the A dependence 
of the one-electron part, and, hence, the charge density 
renormalization due to many-body effects. However, we 
verified that the correction to the total energy due to the 
charge density self-consistency is rather small and within 
our error bars. 

To obtain AF we have computed -^^ — ^^r a discret 
set of values of A ranging from to 1 by performing 
LDA+DMFT simulations with the Coulomb interaction 
scaled accordingly and evaluating {XHij)x in accordance 
with the Migdal formula (we calculated -^^^^^^Ia^o ana- 
lytically, as in this case it is equal to the Hartree-Fock 
approximation to (Hu) computed with the LDA density 
matrix minus the double-counting correction. The result- 
ing value of -^^^I^^Ia^o is small, of order of 0.1 mRy). . 
Then we integrated ^"^^^^^ over A numerically in order 
to obtain AF and its error bar. 

The resulting many-body correction to the electronic 
free energy is displayed in Fig. |3] for temperatures of 
2900 K and 5800 K. There we also show the correspond- 
ing correction to the total energy AE = Edmft^Ei^oa, 
where Ejdmft was computed in accordance with Eq. (3) 
of Ref. hj. Within our error bars the magnitude of AF 
is the same for bcc and hep Fe, which are suggested as 
stable phases of ironQ] and iron-based alloys [H,]!! at the 
Earth's inner core conditions. The magnitude of AF is 
at least several mRy smaller in the case of fee Fe, show- 
ing that the many-body correction may significantly af- 
fect relative energy differences among iron phases at the 
Earth core conditions. One may also notice that the en- 
tropic contribution TAS = AE — AF becomes much 
more significant at the higher temperature, and its con- 
tribution is almost twice larger in the case of the bee 
phase compared with two others. 

The apphcation of the LDA+DMFT theory has the 
most important consequences for the understanding of 
magnetic properties of Fe at the Earth's core conditions. 
In Fig. S] we display the temperature evolution of the 
uniform magnetic susceptibility x the range of tem- 
peratures from 1100 to 5800 K. One may notice that 
in the fee and hep phases the susceptibility exhibits a 
temperature-independent Pauli behaviour expected for a 
FL (a small decrease in x at lower temperatures observed 
in hep Fe is due to a dip in its one-electron DOS in the 
vicinity of the Fermi level, see SI Fig. SI.). In contrast, x 
of bcc Fe features a clear inverse-temperature dependence 
and can be very well described by the Curie- Weiss (CW) 

law X = with /ie// =2.6 /.is and O =1396 K (see 

inset in Fig. |?|). 

One may relate the apparent CW behaviour of the 
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FIG. 3: (Color online) Many-body correction to the total 
(black dashed line) and free (red solid line) energy for the 
three phases of Fe at the volume of 7.05 A/atom at T=5800 K 
(upper panel) and 2900 K (lower panel). The error bars are 
due to the CT-QMC stochastic error. 



magnetic susceptibility in bcc Fe to the high peak at Ep 
present in its LDA DOS (SI Fig. SI), which can lead 
to a strongly temperature-dependent Pauli (band) sus- 
ceptibility . We have computed the Stoner-enchanced 
Pauh susceptibilities Xst = Xo/(l — I * Xo): for all three 
phases, where / is the Stoner parameter and xo is the 
bare uniform Pauli susceptibility computed from the cor- 
responding LDA densities of states. We fitted / to re- 
produce the corresponding values of LDA+DMFT x at 
T =3800 K, the resulting values of / are 0.44, 0.53, and 
0.54 eV in the bcc, fee, and hep phases, respectively. 
Obtained Xst reproduce very well the LDA+DMFT mag- 
neitic susceptibilities of fee and hep Fe, thus confirming 
the FL nature of these phases. In 6cc Fe Xst describes 
weU the CW-hke behaviour of LDA+DMFT susceptibil- 
ity in the range from 3000 to 5800 K. However, Xst devi- 
ates from the LDA+DMFT susceptibility significantly at 
lower temperatures T, which are small compared to the 
characteristic width of the peak at Ep in the LDA DOS, 
see Fig. IH 

An alternative and more interesting source for the ap- 
parent CW behaviour of the uniform susceptibility in bee 
Fe can be a local magnetic moment surviving in this phase 
up to Earth core temperatures. Existence of a local mag- 
netic moment of the constant magnitude also provides a 
natural explanation for the inverse bcc Fe susceptibility 
exhibiting the same linear temperature dependence in the 



4 




3000 4000 
T(K) 



FIG. 4: (Color online) The uniform magnetic susceptibility 
in paramagnetic state versus temperature. The error bars 
are due to the CT-QMC stochastic error. The dashed lines 
with corresponding filled symbols are fits to the enhanced 
Pauli law, see the text. Inset: the inverse uniform magnetic 
suscptibility of hcc Fe is shown in red (empty circles), the 
blue dot-dashed and green (filled circles) lines are fits to the 
Curie- Weiss and enhanced Pauli law, respectively. 



whole range from 1100 to 5800 K, with no significant de- 
viations or noticable peculiarities (see inset in FigU]). 

The value of calculated uniform susceptibility in SI 
units at temperature T=5800 K is equal to 1.7h<10~'*, 
2.0*10"'*, and S.S^IO"** for hep, fee, and hec Fc, respec- 
tively. We would like to underline that our calculated 
uniform magnetic susceptibilities in all three phases of 
Fe are sufficiently high to be important for models of the 
Earth core dynamics and gcodynamo (sec e.g. Rcf. p^). 
An inner core with a paramagnetic susceptibility in the 
range 10"^ - 10"^ SI units, and with a paramagnetic 
relaxation time acting slower than field changes coming 
from the outer core, could also attenuate short frequency 
fiuctuations and become an important factor stabilizing 
the geodvnamojisj. just as an electrically conducting in- 
ner core could stabilize the geodynamo because the inner 
core would have a magnetic diffusion constant indepen- 
dent of the outer core[24l |25|. 

In conclusion, we have carried out a theoretical inves- 
tigation of the role of electronic correlations in the hce, 
fee, and hep phases of Fe at the inner Earth core condi- 
tions using a fully self-consistent LDA-I-DMFT approach. 
We have found that the fee and hep phases remain in 
a Fermi-liquid state, while hcc Fe features a non-Fermi- 
liquid behaviour. We have evaluated a correction to the 
electronic free energy due to many-body effects and found 
that it affects significantly the rellative free energy differ- 
ences, penalizing the fee phase. Most interestingly, our 
results suggests that a local magnetic moment may exist 
in the hcc phase at the inner core conditions and that 
magnetic susceptibilities in all three phases of Fc arc suf- 
ficiently high to stabilize geodynamo. Thus, new models 
of geodynamo as well as the core structure and elasticity 



should include the magnetism of the Earth's core, the 
effect of which so far was never considered. 



METHODS 

In our LDA-l-DMFT calculations Wannier-like functions for 
the Fe-3d shell were constructed by projecting local orbitals 
onto a set of FLAPW Bloch states located within the energy 
window from -10.8 to 4 eV relative to the Fermi level Ef 
(details of the projection procedure can be found in Ref. [ll[). 

We then introduced the calculated local Coulomb interac- 
tion in the density-density form acting between those Wannier 
orbitals. In order to evaluate the strength of the on-site elec- 
tron repulsion on the Fe 3ci shell we employed the constrained 
random-phase-approximation (cRPA) method [2^ . [27j . The 
calculated Coulomb ( U) and exchange ( J) interaction matri- 
ces are well approximated by a spherically-symmetric form 
used in the subsequent calculations, with the parameter U 
(the Slater parameter _Fo) equal to 3.15, 3.04, and 3.37 eV 
and the Hund's rule coupling J equal to 0.9, 0.9 and 0.93 eV 
for bcc, fee, and hep phases, respectively. These results are 
in general agreement with previous calculations of U in com- 
pressed Fe [2^ . We employed the around mean-field form [2^ 
of the double counting correction term throughout. 

The resulting many-body problem has been treated within 
the DMFT framework with the quantum impurity problem 
solved by the exact Continuous-time hybridization expansion 
Quantum Monte-Carlo (CT-QMC) method ^ using 5*10^ 
CT-QMC moves with a measurement performed after each 
200 moves. After completing the DMFT cycle we calculated 
the resulting density matrix in the Bloch states' basis, which 
was then used to recalculate the charge density in the next 
iteration as described in Ref. p^ . To obtain the spectral func- 
tion at the real axis we employed a stochastic version of the 
Maximum Entropy method [3l[ for analytical continuation. 
In order to compute the magnetic susceptibility in uniform 
fields we performed LDA-I-DMFT simulations with the Kohn- 
Sham eigenstates split by the magnetic field H — O.OOSeV/fiB 
directed along the z axis. 
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FIG. SI: Partial LDA densities of states (PDOS) for the irreducible representations of Fe 3d states 
for the three phases at the volume 7.05 A'^/atom. The large peak in the vicinity of Ep in bcc Fe eg 
PDOS is due to a van Hove singuliarity. 
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FIG. S2: Imaginary part of the hybridization function A as function of imaginary frequency for the Cg 
and t2g states at the first DMFT iteration. One may clearly see a decrease in |53A(i(j;)| at a; < 2 eV 
of the Cg states in bcc Fe. However, at higher energies (a; > 5 eV) the eg hybridization function 
decays slowly and becomes larger than the t2g one. The overall one-electron kinetic energies of Cg 
and t2g states in bcc Fe have similar values as explained in the main text. In contrast, in the fee 
phase |A| grows monotonously with decreasing u for both eg and t2g 



